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We describe an experiment to search for a new vector boson A' with weak coupling a' > 6 x 10~*q: 
to electrons (a = e'^/An) in the mass range 65 MeV < tua' < 550 MeV. New vector bosons with 
such small couplings arise naturally from a small kinetic mixing of the "dark photon" A' with the 
photon — one of the very few ways in which new forces can couple to the Standard Model — and 
have received considerable attention as an explanation of various dark matter related anomalies. 
A' bosons are produced by radiation off an electron beam, and could appear as narrow resonances 
with small production cross-section in the trident e"*'e~ spectrum. We summarize the experimental 
approach described in a proposal submitted to Jefferson Laboratory's PAC35, PR-10-009 [T]. This 
experiment, the A' Experiment (APEX), uses the electron beam of the Continuous Electron Beam 
Accelerator Facility at Jefferson Laboratory (CEBAF) at energies of « 1-4 GeV incident on 0.5—10% 
radiation length Tungsten wire mesh targets, and measures the resulting e'^e~ pairs to search for 
the A' using the High Resolution Spectrometer and the septum magnet in Hall A. With a ~ 1 
month run, APEX will achieve very good sensitivity because the statistics of e'^e~ pairs will be 
~ 10, 000 times larger in the explored mass range than any previous search for the A' boson. These 
statistics and the excellent mass resolution of the spectrometers allow sensitivity to a' /a one to 
three orders of magnitude below current limits, in a region of parameter space of great theoretical 
and phenomenological interest. Similar experiments could also be performed at other facilities, such 
as the Mainz Microtron. 



I. INTRODUCTION 

The development of the Standard Model of particle 
interactions is the culmination of a century of searches 
and analyses with fixed-target and colliding beam exper- 
iments. Interactions with new forces beyond the Stan- 
dard Model are currently limited by well-tested gauge 
symmetries to a handful of possibilities. One of the 
few remaining ways for interactions with new sub-GeV 
vector-like forces to arise is for charged particles to ac- 
quire millicharges, ee, under these forces. This occurs 
through a simple and generic mechanism proposed by 
Holdom [2 , in which a new vector particle mixes via 
quantum loops with the Standard Model photon. MeV- 
GeV masses for the A' gauge boson are particularly well- 
motivated in this context. Such sub-GeV forces are a 
common feature of extensions of the Standard Model, 
but existing constraints are surprisingly weak, with lim- 
its at ee < (0.3 - 1) x IQ-^e. 

Fixed-target experiments with high-intensity electron 
beams and existing precision spectrometers are ideally 
suited to explore sub-GeV forces by probing reactions in 
which a new A' vector particle is produced by radiation 
off an electron beam O |4] . The A' can decay to an elec- 
tron and positron pair and appears as a narrow resonance 
of small magnitude in the invariant mass spectrum. The 
production rate of A's, the luminosity, and the mass res- 



olution attainable at, for example, Jefferson Laboratory 
(JLab), the Mainz Microtron, and the SLAG National 
Accelerator Laboratory vastly exceeds what is currently 
available using colliding electron beam facilities. In [3], 
several fixed-target experimental strategies were outlined 
to search for new sub-GeV vector interactions. In this pa- 
per, we summarize a concrete A' search using Jefferson 
Laboratory's Gontinuous Electron Beam Accelerator Fa- 
cility (GEBAF) and the High Resolution Spectrometers 
(HRS) in Hall A [T], highlighting the features that are 
applicable to similar experimental facilities. This exper- 
iment, the A' Experiment (APEX), can probe charged 
particle couplings with new forces as small as 2 x lO^^e 
and masses between 65 MeV and 550 MeV — an im- 
provement by more than two orders of magnitude in cross 
section sensitivity over all previous experiments. 

Fixed-target experiments of this form are particularly 
timely in light of a series of recent anomalies from terres- 
trial, balloon-borne, and satellite experiments that sug- 
gest that dark matter interacts with Standard Model par- 
ticles. Much of this data sharply hints that dark matter 
is directly charged under a new force mediated by an A' 
and not described by the Standard Model. Theoretical 
as well as phenomenological expectations suggest an A' 
mass rriA' ^ 1 GeV and ee < 10~^e. 

In this paper, we shall focus on a search for new vec- 
tor bosons. However, it should be emphasized that this 



experiment will provide a powerful probe for any new 
particle — vector, pseudo- vector, scalar, or pseudo-scalar 
— that has sub-GeV mass and couples to electrons (for 
other collider, accelerator, and direct and indirect as- 
trophysical probes see [5^125): a proposal for an electron 
beam incident on a diffuse Hydrogen gas target using 
the Jefferson Laboratory's Free Electron Laser has been 
discussed in [H]). 



A. Brief overview of the experimental strategy 

The goal of the experiment is to measure the invari- 
ant mass spectrum of electron-positron pairs produced by 
electron scattering on a high-Z target, and search for a 
narrow peak with width corresponding to the instrumen- 
tal resolution. The electron and positron are detected 
in magnetic spectrometers with acceptance over a small 
range of particle momentum and angle, such that each ex- 
perimental setting is sensitive to a mass window ^ ±30% 
about a central mass value. Using four beam energies 
from 1-4 GeV, APEX will scan the e+e" spectrum in 
the mass range 65 MeV to 550 MeV. 

Optimal sensitivity for these masses is achieved by 
studying symmetric e+e^ kinematics, where each particle 
carries approximately half the beam energy and has an 
opening angle « 5 deg relative to the beam. Such small 
effective angles for the spectrometer can be achieved us- 
ing a septum magnet [11 Without a septum magnet, 
lower beam energies and correspondingly wider angles 
could be used to probe the same mass range. The impact 
of the geometry on the physics reach will be reviewed in 



[111 and was discussed in detail in [3]. 



The experimental sensitivity is determined by statistics 
and mass resolution. Given the precision of spectrome- 
ters used, the latter is limited by multiple scattering in 
the target material. In APEX, a long, tilted wire mesh 
target is used to obtain excellent relative mass resolution 
of 0.5%. In addition, different segments of the target will 
enter the spectrometers for different central angles, in- 
creasing the size of mass window probed by each setting. 

With a beam of 80 /^A on 0.5%-10% radiation-length 
targets at various beam energies, we expect to collect true 
coincidence e+e" events with a rate in the range 100-500 
Hz (the expected background and accidental coincidence 
rates within a 2 ns timing window are about an order 
of magnitude lower). The total e^e~ sample size will 
exceed 10® pairs in a 6-day period for each setting, or a 
12-day period for the 4 GeV setting. 

While this paper reflects an experimental setup opti- 
mized for the equipment in Hall A at JLab, many of the 
experimental considerations are also applicable for equip- 
ment available at the Mainz Microtron, JLab Hall B, and 
other experimental facilities. 



B. Expected reach and impact 

APEX will be sensitive to new gauge bosons with cou- 
plings as small as a' /a ^ (6 ~ 8) x 10^® for masses 
in the range 65 — 300 MeV, and couplings as small as 
a' /a - 2 X lO"'^ for larger m^/ < 525 MeV. This is 
about a factor of ~ 3 — 35 times lower in e than exist- 
ing constraints (which assume that the A' couples also 
to muons), and corresponds to ^ 10— 1000 times smaller 
cross-sections. 

The precise mass range probed by this type of experi- 
ment can be varied by changing the spectrometer angu- 
lar settings and/or the beam energies. Thus, other ex- 
perimental facilities may be able to perform experiments 
similar to APEX, but targeting complementary regions 
of parameter space. 

The parameter range probed by APEX is interesting 
for several reasons. This region of mass and coupling is 
compatible with A"s explaining the annual modulation 
signal seen by the dark matter direct detection exper- 
iment DAMA/LIBRA, and also with dark matter an- 
nihilating into yl"s, which explains a myriad of recent 
cosmic-ray and other astrophysical anomalies (see S II B I . 



In addition, and independently of any connection to dark 
matter, the proposed experiment would be the first to 
probe A's of mass > 50 MeV with gauge kinetic mixing 
below e ~ 10"'^, the range most compatible if the Stan- 
dard Model hypercharge gauge force is part of a Grand 
Unified Theory. 

The importance for fundamental physics of discovering 
new forces near the GeV scale cannot be overstated. 



C. The organization of this paper 

The paper is organized as follows. In plj we present 
the physics of hypothetical A' particles, motivation for 
their existence, current limits, and estimated sensitivity 
for potential future analyses of existing data. In §III[ 
we describe A' production in fixed-target experiments. 
In \ IV we describe the experimental setup. In SjV) we 



present the parametrics and the Monte Carlo (MC) simu- 
lations of the QED e^e~ pair production rate and the A' 
signal rate in the proposed setup. We also describe how 
we made the sensitivity plots. Other background rates, 
such as TT"*" or e"*" singles and accidental e^e~ pairs, are 
discussed in |VI[ The expected sensitivity is discussed in 
§VII| The paper is summarized in §VIII| Three appen- 
dices discuss the form factors used to calculate the signal 
and background rates (ijA|), the mass resolution (f|B|, and 
the validation of the rates we obtain with the various MC 
simulations (j|C|). 



II. PHYSICS 

We consider new sub-GeV mass vector bosons — 'dark 
photons' A' — that couple very weakly to electrons 
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(as mentioned previously, similar considerations apply 
to pseudo-vectors, scalars, and pseudo-scalars with sub- 
GeV mass that couple to electrons). It is useful to param- 
eterize the coupling g' of the A' to electrons by a dimen- 
sionless e = g' /e, where e is the electron charge. Cross- 
sections for A' production then scale as a' /a — e^, where 
a' = g'"^ I (47r) and a = / (47r) are the fine-structure con- 
stants for the dark photon and ordinary electromagnetic 
interactions, respectively. This experiment will search 
for A' bosons with mass wia' ~ 65 MeV - 550 MeV and 
ct' /a. > 6 X 10^*, which can be produced by a reaction 
analogous to photon bremsstrahlung (see ^III) and de- 



cays promptly to e^e~ or other charged particle pairs. 
We refer the reader to Figure [l] for a summary of the 
reach of this experiment. 




FIG. 1: Anticipated 2a sensitivity in a' /a = for the A' 
experiment (APEX) at Hall A in JLab (thick blue line), with 
existing constraints on an A' from electron and muon anoma- 
lous magnetic moment measurements, Oe and (see (27)'). 
the BaBar search for T(35') — >■ "ffj,'^jj.~ [28|, and three beam 
dump experiments, E137, E141, and E774 [29H3T] (see [3]). 
The and T(3S) limits assume equal-strength couplings to 
electrons and muons. The red region indicates the region of 
greatest theoretical interest, as described in the text. The 
gray dashed line indicates the scale used for other plots in 
this paper. The irregularity of the reach is an artifact of com- 
bining several different run settings (see Table [ll|. The precise 
mass range probed by this type of experiment can be varied by 
changing the spectrometer angular settings and/or the beam 
energies. We stress this point as other experimental facilities 
may be able to perform experiments similar to APEX, but 
targeting complementary regions of parameter space. 



A. Motivation for New Physics Near the GeV Scale 

New light vector particles, matter states, and their as- 
sociated interactions are ubiquitous in extensions of the 
Standard Model pi However, the symmetries of 

the Standard Model restrict the interaction of ordinary 
matter with such new states. Indeed, most interactions 
consistent with Standard Model gauge symmetries and 
Lorentz invariance have couplings suppressed by a high 
mass scale. One of the few unsuppressed interactions is 
the coupling of charged Standard Model particles tp 



6C = g'A'Jr^ 
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to a new gauge boson A\ which is quite poorly con- 
strained for small g' (see Figure [l]) [5] . Similar couplings 
between the A' and other Standard Model fermions 
are also allowed, with relations between their couplings 
(anomaly cancellation) required for the A' gauge symme- 
try to be quantum-mechanically consistent. For example, 
the A' can couple only to electrons and muons, with op- 
posite charges g'^ — —g'^ ( a [/(l)e-^ boson), or can have 
couplings proportional to the electromagnetic charges qi 
of each fermion, gi = eeqi. 

A' couplings to Standard Model matter with the lat- 
ter structure can be induced by ordinary electromagnetic 
interactions through the kinetic mixing interaction pro- 
posed by Holdom [5], 
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where i^^^ 



d^Al 



dvA'^ is the field strength of the 
A' gauge boson, and similarly Fy'^ is the hypercharge 
field strength. This effect is generic, ensures that the 
A' interactions respect parity, and (as we discuss below) 
naturally produces small g' and A' masses near the GeV 
scale. This mixing is equivalent in low-energy interac- 
tions to assigning a charge eeqi to Standard Model parti- 
cles of electromagnetic charge g^, where e — ey /{cos Ow) 
and 9w is the Weinberg mixing angle. The A' couplings 
to neutrinos and parity-violating couplings are negligible 
compared to Z-mediated effects (see e.g. [IS])- 

As noted in a new gauge boson A' that does not 
couple to Standard Model matter at a classical level 
can still couple through quantum-mechanical corrections. 
For example, loops of any particle X that couples to both 
the A' and Standard Model hypercharge generates mix- 
ing of the form ([2]), with 

[a' /a - 10-^ - 10-4). (3) 



10"^ - 10" 



These quantum effects are significant regardless of the 
mass mx of the particle in question, which could be well 
above the TeV scale (or even at the Planck scale) and 
thus evade detection. 

Smaller e are expected if nature has enhanced sym- 
metry at high energies. For example, it has been con- 
jectured that the strong and electroweak gauge groups 
of the Standard Model are embedded in a grand unified 



3 



theory (GUT) with gauge group SU{5) or larger that is 
broken spontaneously at a high scale Mq ~ 10^^ GeV. 
In this case the mixing ^ is suppressed, 

2 

ecuT 7°^ In (Mg/Mx) ^ IQ-^ - IQ-^, (4) 

where are gauge couplings, e of this size leads to 
effective couplings 

a'/a - 10"* - 10^^ (5) 

As shown in Figure [T] no experiment to date has probed 
the range ofe expected in grand unified theories for tua' ^ 
50 MeV. (From string theory, the possible range of e is 
much larger, - IQ-^^ - IQ-^ [35H3H]-) 

An A' mass near but beneath the weak scale is partic- 
ularly well-motivated, as U{iy symmetry-breaking and 
the resulting A' mass may be determined by the same 
physics that generates the W and Z masses [H]. The 
best candidate for the origin of the weak scale is low- 
energy supersymmetry. In this case, the A' can naturally 
acquire mass suppressed by a loop factor or by ^/e com- 
pared to the weak scale, leading to MeV to GeV-scale 
A' masses [T31 |3S1 HTHii] . In supersymmetric models, 
the gauge kinetic mixing ^ is accompanied by quartic 
interactions 

5£^'j-gYgD\<l>Dm^ (6) 

between the Standard Model Higgs doublet h and any 
scalar c/jd charged under U{1)\ where gy and gu are the 
gauge couplings of Standard Model hypercharge and the 
A' coupling to 4>D, respectively. Electroweak symmetry 
breaking gives h a weak-scale vacuum expectation value, 
so that ^ generates a mass term for (po- For positively 
charged (pn, and sufficiently small bare mass, this mass 
term is negative and triggers J7(l)' breaking by the Higgs 
mechanism. The resulting induced mass for the A' is 

ruA' - yfej^^mw - MeV-GeV, (7) 

V 92 

where 32 is Standard Model SU{2)l gauge coupling and 
ruw is the W-boson mass. The resulting mass is precisely 
in the 50— 1000 MeV range targeted by this experiment. 
Given our e sensitivity, we expect to probe the portion 
of this parameter space with small g^. For example, for 
go ~ 0.04 and e ~ 5 X 10"'' {a' fa ~ 2.5 x 10"^), we 
have ruA' ^ 400 MeV, which can definitively be probed 
by the proposed experiment. Note that the mechanism 
of [/(I)' breaking above does not rely on supersymmetry, 
as any quartic interaction of the form ([6| , with arbitrary 
coupling, can transmit electro- weak masses to the A' . 
Thus, the mass relation ([t]) should not be interpreted 
too literally. 

We stress that the mass of the A' breaks any apparent 
symmetry between it and the photon: though Standard 
Model particles have induced e-suppressed charges under 
the A', any new matter charged under the A' would not 



have any effective coupling to the photon, and would have 
gone undetected. 

An electron beam scattering on a high-Z target such 
as Tungsten will produce A"s through bremsstrahlung 
reactions with a cross-section 




several orders of magnitude larger than in colliding elec- 
tron and hadron beams [7]. The A' can decay to elec- 
trons, and is therefore visible as a narrow resonance in 
the trident e+e^ mass spectrum. 

Such a new gauge boson would constitute the first dis- 
covery of a new gauge force since the observation of Z- 
mediated neutral currents. Besides the obvious phys- 
ical interest of a fifth force, the A' like the Z could 
open up a new "sector" of light, weakly coupled parti- 
cles whose spectrum and properties could be measured 
in fixed-target experiments and flavor factories. The A' 
sector would provide a new laboratory for many physical 
questions, and would be revealing precisely because its 
interactions with Standard Model particles are so weak. 
In particular, if nature is approximately supersymmet- 
ric near the TeV scale, the mass scale of supersymmetry 
breaking for the A' sector is naturally suppressed by e 
times gauge couplings. In this case, supersymmetry could 
be studied easily in the A' sector, and possibly even dis- 
covered there by relatively low-energy experiments before 
Standard Model superpartners are seen at colliders. 

B. Motivation for an A' from Dark Matter 

Dark matter interpretations of recent astrophysical 
and terrestrial anomalies provide an urgent impetus to 
search for A'^s in the mass range 50 MeV - 1 GeV, with 
a coupling e ^ 10^** - 10^^. 

The concordance model of big bang cosmology — the 
"Lambda Cold Dark Matter" (ACDM) model — ex- 
plains all observations of the cosmic microwave back- 
ground, large scale structure formation, and supernovae, 
see e.g. 15 1!) . This model suggests that Standard 
Model particles make up only about 4% of the energy 
density in the Universe, while "dark energy" and "dark 
matter" make up 74% and 22%, respectively, of the Uni- 
verse's energy density. The concordance model does not 
require dark matter to have any new interactions beyond 
gravity with Standard Model particles. However, an in- 
triguing theoretical observation, dubbed the WIMP mir- 
acle, suggests that dark matter does have new interac- 
tions. In particular, if dark matter consists of ^100 GeV 
to 10 TeV particles interacting via the electroweak force 
("weakly interacting massive particles" or "WIMPs"), 
they would automatically have the right relic abundance 
observed today. 

In addition to the WIMP miracle, evidence from 
cosmic-ray data and the terrestrial direct dark matter de- 
tection experiment DAM A/LIBRA strongly suggest that 
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dark matter interacts with ordinary matter not just grav- 
itationally. While the WIMP miracle suggests that dark 
matter is charged under the Standard Model electroweak 
force, we will see that these observations provide impres- 
sive evidence for dark matter interacting with ordinary 
matter through a new force, mediated by a new 50 Me V 
- 1 GeV mass gauge boson. In addition to explaining 
any or all of these observations, dark matter charged un- 
der this new force automatically has the correct thermal 
relic abundance observed today by virtue of its interac- 
tions via the new force carrier, reproducing the success 
of the WIMP dark matter hypothesis. 

The satellites PAMELA [SU] and Fermi [51], the 
balloon-borne detector ATIC [52] , the ground-based tele- 
scope HESS [531 [Mj, as well as other experiments, ob- 
serve an excess in the cosmic-ray flux of electrons and/or 
positrons above backgrounds expected from normal as- 
trophysical processes. If their source is dark matter anni- 
hilation or decay, synchrotron radiation from these elec- 
trons and positrons could also explain the "WMAP haze" 
near the Galactic center [55], which consists of an ex- 
cess seen in the WMAP Cosmic Microwave Background 
data. In addition, starlight near the Galactic center 
would inverse Compton scatter off the high energy elec- 
trons and positrons and produce an excess in gamma- 
rays. A detection of a gamma-ray excess towards the 
Galactic center region in the gamma-ray data obtained 
with the Fermi satellite was recently reported in [56| , and 
has been dubbed the "Fermi haze" . 

Taken together, these observations by several experi- 
mental collaborations provide compelling evidence that 
there is an unexplained excess in cosmic-ray electrons 
and positrons in our Galaxy. Given the firm evidence 
for a 22% dark matter content of the Universe, a very 
natural source of these excesses is dark matter annihila- 
tion. However, two features of these observations are in- 
compatible with annihilation of ordinary thermal WIMP 
dark matter. They instead provide impressive evidence 
that dark matter is charged under a new 17(1)' and anni- 
hilating into the A' , which decays directly into electrons 
and positrons, or into muons that decay into electrons 
and positrons, see Figure [2] (left) (see e.g. [S] [57M5] ). 
These two features are: 

• The annihilation cross-section required to explain 
the signal is 50-1000 times larger than the ther- 
mal freeze-out cross-section for an ordinary WIMP 
that is needed to reproduce the observed dark mat- 
ter relic density. This can be explained if dark 
matter interacts with a new long range force me- 
diated by an ©(GeV) mass gauge boson, which 
allows the dark matter annihilation cross-section 
i{crv)) to be enhanced at low dark matter velocities, 
i.e. (av) (X 1/v. In this case, in the early Universe 
when the dark matter velocity was high (^ 0.3c), 
the annihilation cross-section that determines the 
relic abundance can naturally be the same as that 
of an ordinary WIMP and reproduce the WIMP 
miracle. However, in the Milky Way halo now, the 



dark matter has a much lower velocity {v ~ 10~^c), 
leading to a large increase in the annihilation cross- 
section that is required to explain the cosmic-ray 
data. The enhancement at low velocities through a 
new long-range force is very well known and called 
the Sommerfeld effect j64j . 

• The PAMELA satellite did not see an anti-proton 
excess [5S] , which strongly suggests that dark mat- 
ter annihilation is dominantly producing leptons, 
and not baryons. If dark matter is interacting via a 
0{GeV) mass force particle in order to have a large 
annihilation rate via the Sommerfeld mechanism, 
then annihilations into the force carrier automat- 
ically fail to produce any baryons. Kinematically, 
the force carriers cannot decay into baryons, and 
are instead forced to decay into the lighter charged 
leptons. Thus, annihilation products of dark mat- 
ter are leptonic in this case. 

To explain the additional sources of evidence for a new 
GeV scale force, we briefly summarize the consequence 
for dark matter mass spectra that follow from dark mat- 
ter carrying a charge under a new force. If dark mat- 
ter is charged under a non-Abelian force that acquires 
mass, then radiative effects can split all components of 
the dark matter with size, 6 ~ auAm^VDj where ao 
is the non-Abelian fine structure constant and Amwo 
is the splitting of gauge boson masses [57]. Typically, 
these splittings are Amwn ~ ceD'mwn ~ 1 — lOMeV for 
mwn ~ IGeV [57]. Thus, 6 - 100 keV for ao lO^^. 
These splittings are completely analogous to the split- 
tings that arise between the tt^ and tt" from Standard 
Model SU(2) breaking. If instead a non-Abelian force 
confines at a scale ^ GeV, then a heavy-flavor meson 
can be cosmologically long-lived and thus a dark matter 
candidate |66j . Hyperfine interactions can naturally in- 
duce ~ 100 keV splittings of the dark matter particles in 
this case. We emphasize that the GeV scale force carrier 
particles mediate quantum corrections that generate the 
100 keV and 1-10 MeV splittings of dark matter states 

When mass splittings arise. A' mediated interactions 
of dark matter with ordinary matter as well as dark 
matter self-interactions are dominated by inelastic col- 
lisions [57) . The direct dark matter detection experi- 
ment DAM A/LIBRA as weU as the INTEGRAL tele- 
scope provide intriguing evidence for such interactions. 
The DAMA/Nal W] and DAM A/LIBRA [70^ experi- 
ments have reported an annual modulation signal over 
nearly eleven years of operation with more than 8a sig- 
nificance. Modulation is expected because the Earth's 
velocity with respect to the dark matter halo varies as 
the Earth moves around the sun, and the phase of the 
observed modulation is consistent with this origin. A 
simple hypothesis that explains the spectrum and mag- 
nitude of the signal, and reconciles it with the null results 
of other experiments, is that dark matter-nucleus scatter- 



5 



DM DM* 




FIG. 2: Left: Dark matter annihilation into the dark photon A' , which decays into charged leptons such as electrons and/or 
muons, can explain the cosmic-ray electron and/or positron excesses seen by PAMELA, Fermi, ATIC, HESS, and other exper- 
iments. Right: Dark matter scattering into an excited state off nuclei through A' exchange in direct dark matter detection 
experiments can explain the annual modulation signal observed by DAMA/LIBRA, and the null results of other direct detection 
experiments. 



ing is dominated by an inelastic process, 

XN^X* N, (9) 

in which the dark matter x scatters off a nucleus N into 
an excited state x* with mass splitting S w 100 keV 
[67j . The kinematics of these reactions is also remark- 
ably consistent with all the distinctive properties of the 
nuclear recoil spectrum reported by DAMA/LIBRA. In 
addition, the INTEGRAL telescope [7T] has reported a 
SllkeV photon signal near the galactic center, indicating 
a new source of 1-10 MeV electrons and positrons. This 
excess could be explained by collisions of 0(100 GeV-1 
TeV) mass dark matter into ©(MeV) excited states in 
the galaxy [72] — dark matter excited by scattering de- 
cays back to the ground state by emitting a soft e^e~ 
pair. The SllkeV excess then arises from the subsequent 
annihilation of the produced positrons. 

The existence of an A' may also help explain various 
other particle physics anomalies |27j such as the anoma- 
lous magnetic moment of the muon {(g — 2)^) |73j and 
the HyperCP anomaly [74 . 

While these experimental hints provide an urgent mo- 
tivation to look for an A', it is important to emphasize 
the value of these searches in general. There has never 
been a systematic search for new GeV-scale force car- 
riers that are weakly coupled to Standard Model parti- 
cles. Nothing forbids their existence, and their discovery 
would have profound implications for our understand- 
ing of nature. A relatively simple experiment using the 
facilities available at, for example, Jefferson Laboratory 
and Mainz will probe a large and interesting range of A' 
masses and couplings. 



C. Current Limits on Light L*^(l) Gauge Bosons 

Constraints on new A"s that decay to e~^e~ and the 
search reach of an experiment using the spectrometers of 
Hall A at Jefferson Laboratory are summarized in Fig- 
ure [TJ Shown are constraints from electron and muon 
anomalous magnetic moment measurements, Og and 
|27j . the BaBar search for T{3S) — > 7A' — > jfi'^fi^, and 
three beam dump experiments, E137, E141, and E774 
[3]. The constraints from and the BaBar search as- 
sume that the A' couples to muons — this is the case, for 
example, if it mixes with the photon. If it only couples 
to electrons, then the constraints on a' /a and ttia' in the 
region to which the proposed experiment is sensitive are 
weaker than a' /a < 10~^. 

We refer the reader to O [37] for details on exist- 
ing constraints. Here, we briefly review the constraint 
on e~^e~ — >■ jA' — >■ ^fi'^ii~ derived from the BaBar 
search [7S]. If the A' couples to both electrons and 
muons, this is the most relevant constraint in the re- 
gion probed by the proposed experiment. The analy- 
sis of [7S] was in fact a search for T{3S) decays into a 
pseudoscalar a, T{3S) — >■ 7a — > jfi^fj,^, but can be in- 
terpreted as a limit on A' production because the final 
states are identical. Using Cint ^ 30 fb^^ of data con- 
taining - 122 X 10^ T(3S') events, a 90% C.L. upper 
limit of roughly (1 — 4) x 10~^ on the 7/i+/i~ branch- 
ing fraction was found for ttia' ^ 2m^ — 1 GeV. This 
search would thus be sensitive to about ~ 100 — 500 
events with e'^e~ — >■ jA' ^ij,'^fi~. Requiring that 
a{e+e- 7A') x BR{A' M+M") x Ant < 500, where 
BR{A' IJ,+IJ.~) = l/{2 + R{mA')) for niA' > 2m^ with 

7-, cr(e"^e~ ^ hadrons; E—m a/) ^ i. , i i , . 

^ = a(e+e-^;.+M-;g^m^O ' ^^^^ rcscahug the resultmg 
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constraint to represent a 95% C.L. upper bound, we find 
the constraint depicted in Figure [l] For 
this requires a' /a > 10^^, while the constraint weakens 
at higher masses, especially near the p-resonance. See _ 10 

[57] for a comparison of our sensitivity estimate to those 
previously published. 

We caution that systematic uncertainties in the A' 
limit beyond those quoted in |75j may slightly weaken 
the resulting limit, which should therefore be taken as 
a rough approximation unless further analysis is done. 
First, A' production in B-factories is more forward- 
peaked than the T{3S) decay mode considered in [75] . 
so that the signal acceptance is more uncertain. In ad- 
dition, background distributions in [75| are derived from 
smooth polynomial fits to data collected on the T(4S') 
resonance, which is assumed to contain no signal. This 10 
assumption is not correct for A' production, though the 
resulting systematic effects are expected to be small. 
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D. Sensitivity of Potential Searches using Existing 

Data 

Several past and current experiments have data that 
could be used to significantly improve current limits 
on a' /a, as discussed in [H [51 [27]. Here, we esti- 
mate the potential sensitivity of searches in three chan- 
nels (tt^ — 7> J A' — >• 7e+e~, — > ijA' — ^ yye+e^, and 
e+e~ — )■ J A' — > 7/i+/x~), considering only the statis- 
tical uncertainties and irreducible backgrounds. These 
are likely overestimates, as we are unable to include ei- 
ther systematic uncertainties or significant instrumental 
backgrounds such as photon conversion in the detector 
volume. 

BaBar, BELLE, and KTeV (E799-II) have produced 
and detected large numbers of neutral pions, of order 
10^°, of which roughly 1% decay in the Dalitz mode 
tt" — ;> e^e^7. These experiments can search for the 
decay tt" — > jA' induced by A'-photon kinetic mix- 
ing, which would appear as a narrow resonance over 
the continuum Dalitz decay background. KTeV has the 
largest tt*^ sample, and its e^e^ mass resolution can 
be approximated from the reported measurement of the 
7r° — >■ e+e" branching fraction [75] to be roughly 2 MeV. 
This paper also reports the measured mass distribution 
of Dalitz decays above 70 MeV, from which we estimate 
potential sensitivity to a' /a as small as 5 x 10~^ for 
70 < rn(e+e~) < 100 MeV, as shown by the orange 
shaded region in Figure [3] 

Similarly, KLOE can search for the decay cj) ~^ V^' > 
likewise induced by A' kinetic mixing with the photon, in 
a sample of 10^° (/)'s. An analysis of this data is ongoing 
[77] . We have taken the blue dashed curve in Figure 
[3] from [3], which assumes that mass resolution am is 
dominated by KLOE's 0.4% momentum resolution. We 
have adjusted the contours from [4] to determine a 2a 
contour and enlarged the bin width used to determine 
signal significance from am in [1] to 2.5am,- Above the 



FIG. 3: Anticipated 2a sensitivity in a' /a = for the A' 
experiment (APEX) at Hall A in JLab (thick blue line), com- 
pared with current limits and estimated potential 2cr sensi- 
tivity for A' searches in existing data (dashed lines), assum- 
ing optimal sensitivity as described in the text. From left 
to right: KTeV tt^ — >■ jA' — >■ 'ye'^e~ (orange dashed curve), 
KLOE (f) — > 77 A' — >■ »7e"^e~ (green dashed curve) and Belle 
e'^e~ — ^ yA' — > yfi^fj.~ (gray dashed curve). Existing con- 
straints are as in Figure [T] 



muon threshold, 4> decays are not competitive with B- 
factory continuum production. 

In addition, BaBar and Belle can search for the con- 
tinuum production mode e~^e~ — > jA' — > y^'^ fi'' in their 
full datasets. For example, an analysis of the Belle T(4S') 
data set would increase statistics by a factor of ^ 24 rela- 
tive to the BaBar T{3S) search that we have interpreted 
as a limit above. We have derived the expected sensi- 
tivity (shown as a black dashed line in Figure p| simply 
by scaling the T(35') estimated reach by ^24. These 
searches have not been extended below the muon thresh- 
old because of large conversion backgrounds. 



III. A' PRODUCTION IN FIXED TARGET 
INTERACTIONS 

A' particles are generated in electron collisions on a 
fixed target by a process analogous to ordinary photon 
bremsstrahlung, see Figure [3] This can be reliably es- 
timated in the Weizsacker- Williams approximation (see 
[51 175H50] ). When the incoming electron has energy Eq, 
the differential cross-section to produce an A' of mass 
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FIG. 4: A' production by bremsstrahlung off an incoming 
electron scattering off protons in a target with atomic number 



with energy Ea' 
da 



= xEq is 



dxd cos 9 A' 



[1-X + 



C/2 

x{l 



x)vfi\, {ElxB\) 



f/2 



(10) 



where Z is the atomic number of the target atoms, 
a ~ 1/137, Qa' is the angle in the lab frame between 
the emitted A! and the incoming electron. 



V{x,dA-) = ElxB\, +TO^, 



1 



+ m„x 



(11) 



is the virtuality of the intermediate electron in initial- 
state bremsstrahlung, and x ~ 0.1 — 10 is the Weizsacker- 
Williams effective photon flux, with an overall factor of 
Z^ removed. The form of x and its dependence on the 
A mass, beam energy, and target nucleus are discussed 
in Appendix |Xj The above results are valid for 

32 



me raA' ^ i?o, xd\, <^ 1. 

For niA' ^ m-e, the angular integration gives 

do- _ SZ'^a^e'^x 
dx 



1 



(12) 



(13) 



"A' \ 3(1-2;)^ 

The rate and kinematics of A' radiation differ from mass- 
less bremsstrahlung in several important ways: 

Rate: The total A' production rate is controlled by 

3 2 

Therefore, it is suppressed relative to pho- 

A' 

2 

ton bremsstrahlung by ^ e^-^^^. Additional sup- 
pression from small x occurs for large niA' or small 

Angle: A' emission is dominated at angles 9a' such that 
U{x,9a') ^ 2U{x,0) (beyond this point, wide- 
angle emission falls as 1 /9\, ) . For x near its median 
value, the cutoff emission angle is 



max 



^UlA' 



rrie m 



3/2 ' 



Er, 



e: 



,3/2 



(14) 



which is parametrically smaller than the opening 
angle of the A' decay products, '-^ tuaiIEq. Al- 
though this opening angle is small, the backgrounds 
mimicking the signal (discussed in ^ VI I dominate 
at even smaller angles. 

Energy: A' bremsstrahlung is sharply peaked at a; « 
1, where U{x,Q) is minimized. When an A' is 
produced, it carries nearly the entire beam en- 
ergy — in fact the median value of (1 — a;) is 



The latter two properties are quite important in improv- 



ing signal significance, and are discussed further in \ VI 



Assuming the A' decays into Standard Model particles 
rather than exotics, its boosted lifetime is 



7CT 



0.8cm / Efi V^O^ 
A^cff VlOGeV jl 



^'lOOMeVy 
niA' ) 



(15) 



where we have neglected phase-space corrections, and 
A'cff counts the number of available decay products. If 
the A' couples only to electrons, N^fi = 1. If the A' 
mixes kinetically with the photon, then A'off = 1 for 
•niA' < 2to^ when only A' — >■ e+e~ decays are pos- 
sible, and 2 -I- R{mA') for to^/ > 2m^, where R = 

"^^|r+e"^^M+M°^^'=^rf ^ ED- ^^"^ ''^'^Ses of e and 
•niA' probed by this experiment, the mean decay length 
^0 ^ 250/im is not significant, but the ability to cleanly 
reconstruct vertices displaced forward by a few cm would 
open up sensitivity to considerably lower values of e. 

The total number of A' produced when electrons 
scatter in a target of T ^ 1 radiation lengths is 



N . 



ZW 
' A ml. ^ 



(16) 



where Xq is the radiation length of the target in g/cm^, 
iVo — 6 X 10^^ mole~^ is Avogadro's number, and A is the 
target atomic mass in g/mole. The numerical factor C ~ 
5 is logarithmically dependent on the choice of nucleus 
(at least in the range of masses where the form-factor 
is only slowly varying) and on uia', because, roughly, 
Xo oc (see [3] and [HI])- For a Coulomb of incident 
electrons, the total number of A'^s produced is given by 



10- 



2 /lOOMeV 

iriA' 



(17) 



The spectrometer efficiency can be estimated from 
Monte Carlo simulation of the signal, discussed in SjV] It 
is quite low in APEX, but of course depends on the pre- 
cise spectrometer settings. For example, for niA' = 200 
MeV, Eo = 3.056 GeV, an angular acceptance window 
of = 0.055 - 0.102 rad and |6'y| < 0.047 rad (corre- 
sponding to an HRS central angle of 4.5°) and a momen- 
tum acceptance oi E — 1.452 — 1.573 GeV for both the 
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Electron, P = EO/2 




Electrons 



FIG. 5: The layout of the experimental setup 
details. 



see text for 



positron and one of the electrons, gives a spectrometer 
efficiency of ^ 0.14%. 



IV. EXPERIMENTAL SETUP 

In this section, we describe the experimental setup of 
the APEX experiment in JLab Hall A. Many of these 
features are also readily adaptable to other experimental 
facilities. 

The APEX experiment will measure the invariant mass 
spectrum of e"'"e~ pairs produced by an incident beam 
of electrons on a tungsten target. The experiment uses 
the two high-resolution spectrometers (HRS) [52] avail- 
able in Hall A at JLab (see Table |l] for design specifica- 
tions), together with a septum magnet constructed for 
the PREX experiment , see Figure [s] The physical 
angle of the HRS with respect to the beam line does not 
go below ^12°, but the septum allows smaller angles to 
be probed down to 4° — 5° by bending charged tracks 
outward. The detector package in each HRS available in 
JLab Hall A includes two vertical drift chambers (VDC), 
the single photo-multiplier tube (PMT) trigger scintilla- 
tor counter ("SO counter"), the Gas Cherenkov counter, 
the segmented high-resolution scintilator hodoscope, and 
the double-layer lead-glass shower counter. 

The electron beam has a current of 80 fxA (correspond- 
ing to ^ 7 C on target per day!), and will be incident on 
a solid target located on a target ladder in a standard 
scattering chamber. The target will be made of tungsten 
wires strung together in a horizontal plane orthogonal to 
the beam direction. The target plane will be mounted at 
an angle of about 10 mrad with respect to the horizontal 
plane. The beam will be rastered by ±0.25 mm in the 
horizontal and ±2.5 mm in the vertical direction to avoid 
melting the target. 

The electron will be detected in the the right HRS 
(HRS-R) and the positron will be detected in the left 
HRS (HRS-L). The trigger will be formed by a coinci- 
dence of two signals from the SO counters of the two arms 



Beam 



zig-zag tilted target 
0.01 mm diameter W wires 




Positrons 



FIG. 6: The top view of the tilted target. The beam is 
rastered over an area 0.5x5 mm^ (the latter is in the ver- 
tical direction). The beam intersects the target in four areas 
spread over almost 500 mm. Pair components will be de- 
tected by two HRS spectrometers at a central angle of ±5°. 
Each zig-zag of the target plane is tilted with respect to the 
beam by 0.5° and consists of a plane of parallel wires perpen- 
dicular to the beam. This reduces the multiple scattering of 
the outgoing e'^e~ pair (produced in a prompt A' decay), as 
described in the text. 



and a coincidence of the signal in the SO counters with 
a signal from the Gas Cherenkov counter of the HRS-L 
(positive polarity arm) . A timing window of 20 ns will be 
used for the first coincidence and 40 ns for the second co- 
incidence. The resulting signal will be used as a primary 
trigger of data acquisition (DAQ). An additional logic 
will be arranged with a 100 ns wide coincidence window 
between signals from the SO counters. This second type 
of trigger will be prescaled by a factor 20 for DAQ, and 
is used to evaluate the performance of the primary trig- 
ger. Most of the DAQ rate will come from events with 
a coincident electron and positron within a 20 ns time 
interval. 

Note that since we want to search for a narrow peak 
in the invariant mass spectrum of e"*'e~ pairs, which re- 
quires a high level of statistical precision, it is especially 
important to have a very small level of systematics and a 
smooth invariant mass acceptance. In [1], we show that 
APEX has these properties. 



A. The long tilted target 

The experiment will utilize the standard Hall A scat- 
tering chamber as it is used by the PREX experiment, 
with a target consisting of a 50-cm-long tilted wire mesh 
plane. The concept of the target is presented in Figures [6] 
and [7] The wires comprising each plane are perpendicu- 
lar to the beam-line. The tilt angle of 10 mrad is sufficient 
to ensure stability of the beam-target geometry, and at 
the same time such a tilt angle is 10 times smaller than 
the central angle to the HRS, which results in a reduc- 
tion of the path length traversed by the produced e'*"e~ 
pairs. The wires comprising of each zig-zag plane are 
spaced so that outgoing e~^e~ pairs coming from prompt 
A' decays inside a wire only travel through a single wire 
(for some configurations, the outgoing e'^e~ pair may not 
have to traverse any wire if the A' does not decay inside 
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Target ~ 4-5 zig-zag segments 
zig-zag segment lengtii ~ 10 cm 

Tungsten Wires 

beam spot ~ 0.5 cm K 
beam < ► ^ o o o « ° 

width i "o"q O O O O; O O O O 

50/im o o o ° , -x,, 
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lOO's iJ,m figure not to scale! 



FIG. 7: A schematic close-up view of the target. The figure 
is not to scale. The target consists of 4-5 zig-zag planes, with 
each plane consisting of tungsten wires strung together. Each 
zig-zag plane is ~ 10 cm long, and lies at an angle of ~ 10 
mrad with respect to the beam line. The tungsten wires have 
a radius of ~ lOfim and are spaced at a distance of ~ lOO's 
fim. While each beam electron can traverse up to ~ 10 wires, 
the production and prompt decay of an A' in a wire produces 
e^e~ pairs that have an angle of ~ m^' /-Ebeam, large enough 
for them to miss the next wire — this greatly reduces the 
multiple scattering, and is the reason for not using a target 
foil. The beam width is ~ 50/^m, which translates into a 
~ 0.5cm large beam spot along the target plane. The vertical 
rastering of the beam of ~ 0.5mm moves the beam spot ~ 5cm 
back-and-forth along the target plane — this helps to prevent 
the beam from melting the target. 

a wire). For wire thickness of ~ 10"'^ radiation lengths, 
this considerably reduces the multiple scattering in the 
target versus that in a true foil and leads to a much bet- 
ter mass resolution. The maximum number of wires that 
a beam electron can pass through per plane is ~ 10 in 
the configuration illustrated in Figure [6] or Figure [7] as- 
suming 10/im diameter tungsten wires. Wires as thick 
as IbfiTu can be used without significantly compromising 
mass resolution. 

The plane of the target wire mesh will be vertical. The 
mesh plane will have 4-5 zig-zags, each with length 5 — 10 
cm, which result in multiple intersection points allowing 
an extra factor of 4-5 for rejection of accidental tracks 
in offline analysis. (Figure |6] illustrates a target with 4 
zig-zags.) 

The central angle of the spectrometer varies with the 
position of the target. In this experiment, such a varia- 
tion is very useful because it extends the range of invari- 
ant mass covered with one setting of the spectrometers. 
For several settings suggested in our run plan, only two 
planes of wire mesh are needed, one at the front and one 
at the back of the acceptance region. 

There are two considerations to take into account when 
selecting the material. The first consideration is to 
achieve the highest possible ratio of signal to background, 
while keeping the background rate low enough so as not 



TABLE I: Main design characteristics of the Hall A High Res- 
olution Spectrometers at nominal target position (see 82, for 
more details). The resolution values are for the full- width at 
half-maximum. These parameters correspond to a point tar- 
get and do not include the effects of multiple scattering in 
the target and windows. In the calculation of the invariant 
mass resolution the effect of multiple scattering in the target 
was taken into account. The vacuum coupling of the scatter- 
ing chamber and the spectrometer allows one to avoid using 
windows. 



to overwhelm the triggering and DAQ system. The sec- 
ond is whether a thin foil or a thin wire of a particular 
material is available. Large backgrounds come from pions 
produced in photo-production from nucleons, and from 
electrons produced in the radiative tail of electron-proton 
elastic scattering. These backgrounds do not mimic the 
signal, but if their rate is too large, they can overwhelm 
the DAQ system. These considerations favor the use of a 
tungsten target, with a total thickness between 0.5% and 
10% radiation length, with thicker targets used in higher- 
energy runs. Reduction of the thickness at low energies 
is required to limit the rate in the electron spectrometer 
and also minimizes the multiple-scattering contribution 
to the pair mass resolution. 

The heat load of the target is also an important con- 
sideration. This is mitigated by rastering the beam and 
using materials like tantalum or tungsten. The tilted tar- 
get cools from a large area; for example, with the raster 
size 0.5 X 5 mm and proposed geometry (Figure [6]), the 
cooling area is 20 cm^. For the parameters of APEX (an 
electron beam of 80 /iA on a 10% Xq tungsten target) 
the head load is about 140 W (or 7 W/cm^), which re- 
sults in the equilibrium target temperature of 1000°K. 
Experimental study has demonstrated that 1 kW/cm^ is 



Configuration 


QQD„Q Vertical bend 


Bending angle 


45° 


Optical length 


23.4 m 


Momentum range 


0.3 - 4.0 GeV/c 


Momentum acceptance 


-4.5% < 5p/p <-(-4.5% 


Momentum resolution 


1x10"'' 


Dispersion at the focus (D) 


12.4 m 


Radial linear magnification (M) 


-2.5 


D/M 


5.0 


Angular range HRS-L 


12.5° - 150° 


HRS-R 


12.5° - 130° 


Angular acceptance: Horizontal 


±30 mrad 


Vertical 


±60 mrad 


Angular resolution ; Horizontal 


0.5 mrad 


Vertical 


1.0 mrad 


Sohd angle at (5p/p = 0, yo = 


6 msr 


Transverse length acceptance 


±5 cm 


Transverse position resolution 


1 mm 
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(b) 



FIG. 8: Sample diagrams of (a) radiative trident (7*) and 
(b) Bethe-Heitler trident reactions that comprise the primary 
QED background to A' — >■ £'^£~ search channels. 



a safe level for a tungsten foil target [83], so we expect 
that the wire mesh target will perform quite well. 



V. SIGNAL AND TRIDENT KINEMATICS 

The stark kinematic differences between QED trident 
backgrounds and the A' signal are the primary consid- 
erations in determining the momentum settings of the 
spectrometers. As we will show in §VI[ QED tridents 
dominate the final event sample after offline rejection of 
accidentals, so we consider their properties in some detail 
here. 

The irreducible background rates are given by the di- 
agrams shown in Figure [Sj These trident events can 
be usefully separated into "radiative" diagrams (Figure 
^a)), and "Bethe-Heitler" diagrams (Figure ^h)), that 
are separately gauge-invariant. 

We have simulated the production of these continuum 
trident background events in QED using the nuclear elas- 
tic and inelastic form- factors in |78j . The simulation is 
done using MadGraph and MadEvent [53] to compute 
the matrix elements for Z (e+e^) Z exactly, 

but neglecting the effect of nuclear excitations on the 
kinematics in inelastic processes. The MadEvent code 
was modified to properly account for the masses of the 
incoming nucleus and electron in event kinematics, and 
the nucleus is assumed to couple with a form-factor G2 
defined in Appendix A. 

The continuum trident background was simulated in- 
cluding the full interference effects between the diagrams 
in Figure [8] In addition, a "reduced-interference" ap- 
proximation simplifies the analysis and is much less com- 
putationally intensive. In this approximation, wc treat 
the recoiling e~ and the e~ from the produced pair as dis- 
tinguishable. Furthermore, we separate trident processes 
into the radiative diagrams (Figure [sj^a)) and the Bethe- 
Heitler diagrams (Figure [8](b)), and we calculate the 
cross-section for both of these diagrams separately. This 
approximation under-estimates the background rates by 
a factor of about 2-3 in the range of A' masses and beam 
energies considered in this paper and [1]. For the reach 
analysis discussed below, we have used differential distri- 
butions computed in the "reduced-interference" approx- 



imation, then rescaled to the cross-section for the full- 
interference process. 

The contribution from the radiative diagrams (Figure 
l8](a)) alone is also useful as a guide to the behavior of 
a' signals at various masses. Indeed, the kinematics of 
the A' signal events is identical to the distribution of 
radiative trident events restricted in an invariant mass 
window near the A' mass. Moreover, the rate of the 
A' signal is simply related to the radiative trident cross- 
section within the spectrometer acceptance and a mass 
window of width 5m by l3] 



da{e- Z ^ e- Z{A' ^ l+i-)) ( iire^ 



( 



ruA' 



da{e-Z e-Z{j* i+£-)) \2N^ifa J \ dm 

(18) 

where N^s counts the number of available decay products 
and is defined below equation (15 1. This exact analytic 



formula was also checked with a MC simulation of both 
the A' signal and the radiative tridents background re- 
stricted to a small mass window Sm, and we find nearly 
perfect agreement. Thus, the radiative subsample can be 
used to analyze the signal, which simplifies the analysis 
considerably. 

It is instructive to compare kinematic features of the 
radiative and Bethe-Heitler distributions, as the most 
sensitive experiment maximizes acceptance of radiative 
events and rejection of Bethe-Heitler tridents. Although 
the Bethe-Heitler process has a much larger total cross- 
sections than either the signal or the radiative trident 
background, it can be significantly reduced by exploit- 
ing its very different kinematics. In particular, the A' 
carries most of the beam energy (see discussion in ^ HI I , 



while the recoiling electron is very soft and scatters to a 
wide angle. In contrast, the Bethe-Heitler process is not 
enhanced at high pair energies. Moreover, Bethe-Heitler 
processes have a forward singularity that strongly favors 
asymmetric configurations with one energetic, forward 
electron or positron and the other constituent of the pair 
much softer. 

These properties are discussed further in the Appendix 
of [3] , and illustrated in Figure |9j which shows a scatter- 
plot of the energy of the positron and the higher-energy 
electron for the signal (red crosses) and Bethe-Heitler 
background (black dots). The signal electron-positron 
pairs are clearly concentrated near the kinematic limit, 
E{e~^) + E{e~) « Ebeam- Background rejection is opti- 
mized in symmetric configurations with equal angles for 
the two spectrometers and momentum acceptance of each 
spectrometer close to half the beam energy (blue box). 

While the signal over background (S/B) can be sig- 
nificantly improved with a judicious choice of kinematic 
cuts, the final S/B in a small resolution limited mass 
window is still very low, ^ 1%. A "bump-hunt" for a 
small signal peak over the continuous background needs 
to be performed. This requires an excellent mass reso- 
lution, which has an important impact on target design 
and calls for a target that is tilted with respect to the 
beam line (see Appendix B for a discussion of the mass 
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Background vs. Signal Kinematics 
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FIG. 9: Positron and electron momenta in A' signal events 
with rriA' = 200 MeV (red crosses) and in Bethe-Heitler back- 
ground events, for a 3 GeV beam energy. Comparably sized 
signal and Bethe-Heitler samples were used to highlight the 
kinematics of both; in fact the expected signals are much lower 
than the Bethe-Heitler process (see Figure 111. The clustering 
of A' events at high momenta near the kinematic limit and 
of Bethe-Heitler events along both axes are evident. A spec- 
trometer acceptance window that optimizes signal sensitivity 
is indicated by the blue box. 



small statistics sample is generated, sufficient to 
obtain the cross-section from MadEvent. 

• The cross-section ratio of the full continuum back- 
ground (with interference effects) to the sum of the 
Bethe-Heitler and radiative tridents is calculated, 
and represents a multiplicative factor by vifhich the 
latter must be multiplied to get the background 
cross-section. 

• The rates of all reactions impinging the spectrome- 
ter acceptance were calculated by integrating over a 
chosen target profile, which usually extended from 
4.5 to 5.5 degrees. For Bethe-Heitler, radiative tri- 
dents, and the continuum trident background, the 
calculation of the rate was performed as a function 



of the invariant mass of the 



pairs. 



• Using the expressions in Appendix B, we calculated 
the mass resolution 5m- We then tiled the accep- 
tance region with bins of size 2.5 x (5™ in invariant 
mass. 

• As a function of a' /a, the total number of signal 
(5) and background (B) events was calculated with 
the help of ( 18 ) for each bin. 



• We then set S/V B = 2, and solved for a' /a. 
This procedure was used to calculate the reach in the 



a'/a and tua' parameter space shown in { VH Further 
improvements may be obtained by more sophisticated 
analysis cuts such as the use of matrix element methods 
(see e.g. [12]). 



resolution). 



VI. BACKGROUNDS 



A. Calculation of the e reach 

For all cross sections and rates of reactions described in 
this paper and J, , Monte Carlo based calculations were 
performed over a grid of beam energy settings and central 
spectrometer angular settings. Interpolation was used to 
extend this grid continuously to intermediate beam ener- 
gies and angles — all rates exhibited expected power law 
behavior, thereby providing confidence in the reliability 
of an interpolation. Additional cross checks at specific 
points were performed to test the accuracy of our inter- 
polation, which was generally better than ^ 5%. 

In order to calculate the a'/a reach for a particular 
choice of target nucleus, spectrometer angular setting, 
profile of wire mesh target, and momentum bite, the fol- 
lowing procedure is performed: 

• Monte Carlo events are simulated for the Bethe- 
Heitler, radiative tridents, and the continuum tri- 
dent background including the full interference ef- 
fects between the diagrams. The latter back- 
ground is computationally intensive, and only a 



In this section, we present the results of an analysis 
of the expected backgrounds for the A' search. Table [III 
summarizes the expected singles rates, trigger rates, and 
coincidence rates. For more details on how we calculated 
the background rates we refer the reader to |lj . 

Important backgrounds come from electron, pion, and 
positron singles. There are three contributions to the 
electron singles rate in the HRS at the proposed mo- 
mentum settings, namely inelastic scattering, radiative 
elastic electron-nuclei scattering, and radiative quasi- 
elastic electron-nucleon scattering. Our calculations of 
the electron, pion, and positron singles rates were checked 
against measurements made by experiment E03-012 for 
a 5 GeV electron beam incident on a hydrogen target, at 
6° 2-GeV HRS setting. The final values of the electron 
and pion rates were obtained by means of the "Wiser" 
code |85j : positron singles rates from trident reactions 
were calculated using MadCraph and MadEvent [53] , de- 
scribed in S|V) 

Using our calculations of the singles rates, we compute 
the rate of accidental coincident triggers arising from an 
6+ in the HRS-R and an e~ in the HRS-L within the 
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Settings 


A 


B 


C 


D 


Beam energy (GeV) 


2.302 


4.482 


1.1 


3.3 


Central angle 


5.0° 


5.5° 


5.0° 


5.0° 


Effective angles 


(4.5,5.5) 


(5.25,6.0) 


(4.5,5.5) 


(4.5,5.5) 


Target T/Xq (ratio") 


4.25% (1:1) 


10% (1:1) 


0.58% (1:3) 10% (1:1) 


Beam current (fiA) 


80 


80 


80 


80 


Central momentum (GeV) 


1.145 


2.230 


0.545 


1.634 


Singles (negative polarity) 










e- (MHz) 


4.5 


0.7 


6. 


2.9 


TV- (kHz) 


640. 


2200 


36. 


2500. 


Singles (positive polarity) 










7r++p (kHz) 


640. 


2200 


36. 


2500. 


e+ (kHz) 


31. 


3.6 


24. 


23. 


Trigger/DAQ: 










Trigger*" (kHz) 


4. 


0.4 


3.2 


3.4 


Signal to background: 










Trident (Hz) 


610 


70 


350 


530 


Two-step (Hz) 


35 


15 


5 


75 


Background'^ (Hz) 


70 


1.3 


70 


35 



" The listed total target thickness is split between two sets of wire mesh planes, located at different 
z to produce the two indicated effective angles. The numbers in parentheses denote the ratio of 
target thickness at the larger effective angle to that at the smaller effective angle. 

Trigger: Coincidence with 20 ns time window between SO-N (assuming pions are rejected by a factor 
of 100) and SO-P signals. 

Dominated by e^e~ accidental rate. We assume pion rejection by a factor of 10* in offline cuts, a 2 
ns time window and additional factor of 4 rejection of accidentals from the target vertex. Further 
rejection using kinematics is expected, but not included in the table. 



TABLE II: Expected counting rates for the proposed experiment. 



trigger timing windows. These accidental coincidences 
are a dominant part of the recorded events in APEX, 
and determine the maximum rate at which potential sig- 
nal trident events can be recorded. A typical composi- 
tion of the single rate in the spectrometers is expected 
to be e~ /7T~ w 80/20 in the negative polarity arm and 
TT~^/p/e^ ~ 80/19/1 in the positive polarity arm. The 
fraction of the true coincidence events could be up to 
50% for the e~7r'^ rate within a 2 ns time window, and 
could be significant for the e~p events in certain regions 
of momenta. 

Besides the trident events discussed in SjV) an addi- 
tional source of true coincidence events is the "two-step" 
(incoherent) trident process, in which an electron radi- 
ates a real, hard photon in the target that subsequently 
converts to a high-mass e~^e~ pair. For thin targets, this 
process is suppressed compared to the trident rate, and 
so it is sub-dominant for all the settings we consider. 

The consideration of these rates determine trigger rates 
and upper bounds on offline accidental rates shown in 
Table nil 



VII. MEASUREMENTS AND REACH IN APEX 

We consider a twelve-day run in the configuration "B" 
of Table |ll] and six-day runs in each of the remaining 
configurations, to search for new resonances in e'^e~ tri- 
dent spectra from 65 to 550 MeV. For settings "A" and 
"C", the target thickness and beam current have been 
optimized to accumulate the largest possible sample of 
trident events without saturating the data acquisition 
system. Settings "B" and "D" are far from data acqui- 
sition limits, but we do not use T/Xq > 10% to avoid 
limits on the total radiation produced (this can possibly 
be side-stepped at other facilities). 

The mass range from 65 to 550 MeV is chosen to take 
advantage of the Hall A HRS spectrometers, as well as 
for its theoretical interest. Lower masses are more ef- 
fectively probed by using lower beam energies, improved 
forward acceptance, and/or vertexing (see also [T^]). Set- 
tings at higher masses are possible but have significantly 
reduced sensitivity and are better suited to exploration 
with higher-acceptance equipment and an experiment op- 
timized to accept muon and pion pairs as well as elec- 
trons. 
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Sensitivity of Proposed Run Plan 
0.1 0.3 0.5 




0.1 0.3 0.5 

e^e'iA') Mass (GeV) 

FIG. 10: Anticipated 2a sensitivity in a' /a = for APEX 
[1] for the settings given in Table [ll] (assuming a six-day run 
in configuration "A" , "C" , and "D" and a twelve-day run in 
"B"). Existing constraints are shown in the gray shaded re- 
gions. The colored curves correspond to the sensitivity in each 
of the individual energy settings, and the thick gray curve re- 
flects the sensitivity of a combined analysis. 



In each setting, the proposed experiment will accumu- 
late between 70 and 300 million trident events. With 
these statistics, it will be possible to search offline for 
small resonances comprising a few thousandths of the 
collected data in a resolution-limited window. This will 
allow sensitivity to new gauge boson couplings a'/a as 
low as 10~^ over the broad mass range probed by APEX, 
as summarized in Figure |10[ This sensitivity would im- 
prove on the cross-section limits from past experiments 
by a factor of - 10 - 1000. 

As a specific example, we have illustrated the expected 
sensitivity of setting "A" to A' signals with different e 
Each component of the target populates 



in Figure 11 



a different invariant mass distribution; for simplicity we 
consider only the contribution from the front planes of 
the target, with 9f,ff « 5.5° (recall that the target is 
extended along the beam line and consists of 4-5 planes 
in a zig-zag configuration). The top panel illustrates the 
absolute size of A' signals at TOy^/ = 200 MeV compared 
to the continuum trident background (gray line) and the 
size of 2 and 5-sigma statistical fluctuations (blue and 
green dashed lines), while the bottom panel illustrates 
how the same signals would appear after subtracting a 
smooth parameterization of the background. The purple 
curves in each panel corresponds to an A' signal with 
a' /a = 7 X 10"^ at 200 MeV, which according to the 
estimates in §11 D| would not be seen or excluded at 2a 
by a future KLOE search in — r/A' . The red curve 
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FIG. 11: Comparison of signal rates in six days of running at 
setting "A" to expected background and statistical sensitivity. 
Top: The resonances in purple and red lines correspond to A' 
signals at 200 MeV, smeared by a Gaussian to model detec- 
tor resolution and multiple scattering, with a' /a = 6.5 x 10~® 
and 1.3 x 10~^, respectively. The upper (purple) signal is just 
beyond the 2<j expected sensitivity of a KLOE analysis, while 
the lower (red) signal corresponds to the "5(t" sensitivity (not 
including a trials factor) of this experiment. The gray line is 
the simulated invariant mass distribution for the continuum 
trident background, and the blue and green dashed lines re- 
flect the size of 2 and 5a Poisson fluctuations. Bottom: The 
gray line corresponds to the bin-by-bin differences between 
pseudodata containing no signal and a smooth fit to this pseu- 
dodata. Analogous subtractions when a signal is present are 
shown in purple and red, with the same a' /a as in the top 
figure. Again the blue and green dashed lines reflect the size 
of 2 and 5a Poisson fluctuations. 



has a' /a = 1.3 x 10 ^, corresponding to the expected 
"5cr" sensitivity (not accounting for the trials factor) in 
APEX. 



VIII. CONCLUSIONS 

This paper summarized a new experiment ( "A' experi- 
ment^ , or "APEX" ) that has been proposed to the Jeffer- 
son Laboratory's PAC 35 [T]. The experiment proposes 
to use 30 days of beam to measure the electron-positron 
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pair mass spectrum and search for new gauge bosons A' 
in the mass range 65 MeV < niA' < 550 MeV that have 
weak couphng to the electron. Parametrizing this cou- 
phng by the ratio a' /a that controls the A' production 
cross-section, this experiment would probe a' /a as small 
as (6 - 8) X 10~^ at masses from 65 to 300 MeV, 
and a' /a ^ (2 — 3) x lO"*" at masses up to 525 MeV, 
making it sensitive to production rates 10-1000 times 
lower than the best current limits set by measurements 
of the anomalous muon magnetic moment and by direct 
searches at BaBar. The experiment uses the JLab elec- 
tron beam in Hall A at energies of about 1, 2, 3, and 4 
GeV incident on a long (50 cm) thin tilted tungsten wire 
mesh target, and both arms of the High Resolution Spec- 
trometer at angles between 5.0° and 5.5° relative to the 
nominal target position. The experiment can determine 
the mass of an A' to an accuracy of ~ 1-2 MeV. 

While this paper was motivated by a specific exper- 
imental proposal for JLab Hall A, very similar experi- 
ments are possible at other experimental facilities, such 
as the Mainz Microtron or JLab Hall B. Many of the 
considerations discussed in this paper are applicable to 
these other facilities. 

Constraints on new vector bosons with mass near 50 
MeV - 1 GeV are remarkably weak. However, such light 
force carriers are well motivated theoretically, and sev- 
eral recent anomalies from terrestrial and satellite exper- 
iments suggest that dark matter interacting with Stan- 
dard Model particles has interactions with new vector 
bosons in precisely this mass range. The proposed ex- 
periment can probe these hypothetical particles with a 
sensitivity that is un-rivaled by any existing or planned 
experiment. 



Acknowledgements 

We thank Nima Arkani-Hamed, James Bjorken, Dou- 
glas Finkbeiner, Mathew Graham, John Jaros, Takashi 
Maruyama, Ken Moffeit, Richard Partridge, Michael Pe- 
skin, Maxim Pospelov, Allen Odian, Art Snyder, Jay 
Wacker, and Neal Weiner for many useful discussions 
and encouragement. We also thank Johan Alwall for an- 
swering questions about Madgraph. We especially thank 
Peter Bosted, Kees de Jager, Doug Higinbotham, John 
LeRose, and Yi Qiang for their help in preparing the pro- 
posal for JLab Hall A. RE and PS are supported by the 
US DOE under contract number DE-AC02-76SF00515. 
RE, PS, and NT would hke to thank the KavH Institute 
for Theoretical Physics in Santa Barbara for hospital- 
ity during part of this research. PS and NT also thank 
the Perimeter Institute for Theoretical Physics for hos- 
pitality. We also thank the SLAG National Accelerator 
Laboratory for providing funding for the "Dark Forces" 
workshop in September 2009. 



Appendix A: Effective Photon Flux, Target Nucleus 
and Beam-Energy Dependence 

In this appendix we summarize the formulas used in 
Section |III| for the reduced effective photon flux x, and 
highlight its dependence on the A' mass, target nucleus, 
and beam energy. The effective photon flux x is obtained 
as in [751 ES] by integrating electromagnetic form- factors 
over allowed photon virtualities: 

For a general electric form factor G2{t), 



X 



dt 



t-tr 



(Al) 



(the other form factor, Gi{t), contributes only a negligi- 
ble amount in all cases of interest). Since we are domi- 
nated by a coherent scattering with G2 oc Z^, it is useful 
to define a reduced photon flux. 



X = x/z' 



(A2) 



and so is logarit 
and tm.ax — ni A, 



The integral in (|A1|) receives equal contributions at all t, 

^2 



imically sensitive to t^. 



W2So) 



For most energies in question, G2{t) is dominated by 
an elastic component 



G2Mt) 



1 + aHj \l + t/d 



1 



(A3) 



where the first term parametrizes electron screening (the 
elastic atomic form factor) with a = 111 Z^^^"^ /nie, and 
the second finite nuclear size (the elastic nuclear form fac- 
tor) with d = 0.164 GeV^A-2/3. We have multiphed to- 
gether the simple parametrizations used for each in |78j . 



The logarithm from integrating ( Al I is large for tmin < d, 



which is true for most of the range of interest. However, 
for heavy A' , the elastic contribution is suppressed and 
is comparable to an inelastic term. 



G2,in{t) — 



1 + a'H 



1 



4m- 



-{A 



(1 



Z, (A4) 



0.71 GcV^ 



where the first term parametrizes the inelastic atomic 
form factor and the second the inelastic nuclear form 
factor, and where a' — 773Z^^/^/me, mp is the proton 
mass, and (Xp = 2.79 [7S]. This expression is valid when 
t/Airip is small, which is the case for uia' in the range 
of interest in this paper. At large t the form factors will 
deviate from these simple parameterizations but can be 
measured from data. One can show that the contribu- 
tion from the other inelastic nuclear form factor Gi{t) is 
negligible. 

The resulting reduced form factor x{m'^,Eo) = x/Z^ 
are plotted in the left panel of Figure [12] as a function 
of e+e" mass for various electron energies (1, 2, 3, and 
4 GeV) incident on a Tungsten target. The relative ef- 
ficiency of A' production in targets of different composi- 
tions but the same thickness in radiation lengths is given 
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FIG. 12: Top; The factor x = x/Z'^ defined in (|A2| and 
(All as a function of e+e~ mass for (bottom to top) 1, 2, 3, 

Bottom: 



and 4 GeV i ncide nt electrons on a Tungsten target. 
The ratio of ( A5 1 A' production rates per radiation length for 
Silicon and Tungsten targets, as a function of invariant mass 
and for beam energies (top to bottom at 0.4 GeV) 1, 2, 3, and 
4 GeV incident electrons. 



by the ratio 

R{Zi,Z2) 



Xo{Zi)x{Zi,t)/AiZ, 

Xo{Z2)x{Z2,t)/A{Z2 



(A5) 



For example the ratio R{Si, W) is shown in the right 
panel of Figure [l2j again as a function of e+e~ mass for 
beam energies between 1 and 4 GeV. 



where 5g is the angular resolution of the electron or 
positron, and Sp/p is the momentum resolution of the 
HRS, which is always less than 3 x 10^'' (in our esti- 
mates for the reach of e, we take Sp/p to be equal to this 
upper bound). We have 



(B2) 



where Shbs is the HRS angular resolution, which is ~ 0.5 
mrad in the horizontal direction and ^ 1 mrad in the ver- 
tical direction. Moreover, S™^ represents the degradation 
of the resolution due to multiple Coulomb scattering in 
the target. It is given by the standard formula ^81) 



13.6 

p[MeY] 



1 + 0.038 In 



(B3) 



where t is the thickness in radiation lengths of the ma- 
terial along the path of the particle, Xq is the radiation 
length of the target in g/cm^, and p is the momentum of 
particle in MeV. 

For the proposed experiment, the thickness of the tar- 
get along the direction of the beam line varies from 
t — 0.003X0 to t = 0.09X9. However, in the case of a 
foil target or a target composed of several thin wires, the 
distance traversed by trident electron-positron pairs can 
be significantly smaller because the electron and positron 
have relatively large angles with respect to the beam line. 
For a foil, we can take t ^tf/ sm{9p — Of), where tf is 
the foil thickness and 6f its angle relative to the beam 
line. In this case the effective thickness traversed by the 
beam is Tq — tf / sin 9 f . 

In the case of a target composed of multiple wires, 
as was assumed in determining the experimental sensi- 
tivity, the wires can be spaced widely enough that the 
pair-produced particles need only travel through a sin- 
gle wire. In this case, t is typically the radius of the 
wire, which for Tungsten wire targets can be as small as 
10/im, or 3 X 10~^Xo. In this case, we find that the HRS 
angular resolution, Shrs: is comparable to the multiple 
scattering (5™* in the proposed experiment. 



Appendix C: Monte Carlo validation with E04-012 
data 



Appendix B: Mass resolution 

In this appendix, we briefly describe an estimate of 
the mass resolution of the spectrometer. Since we are 
looking for a small bump on the invariant mass spectrum 
distribution, an excellent mass resolution is essential to 
obtain a good reach in e. 

The mass resolution of the spectrometer, Sm, is roughly 
given by 



In this appendix, we briefly describe a validation of the 
Monte Carlo (MC) simulation of the signal, Bethe-Heitler 
and radiative trident backgrounds (shown in Figure|8]and 
discussed in j|V|, and the positron singles. 

We first discuss a comparison of the MC with previous 
experimental results from the JLab experiment E04-012 
|86j . This experiment consisted of a 5.01 GeV, lA.BfiA 
electron beam incident on a 1.72% radiation length liquid 
Hydrogen target. The e+ singles rate was measured to be 
~ 1.1 kHz in a momentum window of ±4% around 1.93 
GeV and an angular acceptance of 4.5 msr with an aspect 
ratio of 2-to-l centered at an angle of 6°. The e'^e~ coin- 
cidence rate was measured at ^ 4 Hz for the same angular 
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acceptance for both the electron and positron arm, and 
with a momentum window of ±4% around 1.93 GeV for 
the positron and ±4% around 1.98 GeV for the electron. 
We simulated this with MadGraph and MadEvent [53^ as 
described in using a form factor for Hydrogen given 
in [78]. We find a e"*" singles rate of ^ 965 Hz and an 
e+e" coincidence rate of 3.9 Hz, which agrees with the 
measured rates to within ~ 19% and a few percent, re- 



spectively. 

We have also verified the implementation of form fac- 
tors in Monte Carlo by simulating photo-production of 
electrons and muons off Tungsten and Beryllium with 
MadGraph and MadEvent. The resulting cross-sections 
agree to within 30% with published computations in the 
Weizsacker- Williams approximation (791. 
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